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Nonperturbative renormalisation of composite operators with overlap 
quarks* 
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We compute non-perturbatively the renormalisation constants of composite operators on a 16 3 x 32 lattice 
with lattice spacing a = 0.093 fm for the overlap fermion action by using the regularisation independent (RI) 
scheme. The quenched gauge configurations are generated by tadpole improved plaquette plus rectangle action. 
We test the perturbative continuum relation Za = Zv and Zs ~ Zp and find that they agree well above fi = 
1.6 GeV. We also perform a Renormalisation Group analysis at the next-to-next-to-leading order and convert the 
renormalisation constants to the MS scheme. 
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1. INTRODUCTION 

Lattice QCD is a unique tool with which we can 
compute physical observables non-perturbatively 
from first principles. Renormalisation of lattice 
operators is an essential ingredient needed to de- 
duce physical results from numerical simulations. 
In this paper we study the renormalisation prop- 
erties of composite bilinear operators with the 
overlap quark action. 

In principle, the renormalisation of a quark 
bilinear can be computed by lattice perturba- 
tion theory. However, lattice perturbation theory 
converges slowly and the higher-order corrections 
may not be small, thus introducing a large un- 
certainty in the calculation of the renormalised 
matrix elements in some continuum scheme. 
To overcome these difficulties, Martinelli et al. 
PP have proposed a promising non-perturbative 
renormalisation procedure. The procedure allows 
a full non-perturbative computation of the ma- 
trix elements of composite operators in the Reg- 
ularisation Independent (RI) scheme |1I2| . The 
matching between the RI scheme and MS, which 
is intrinsically perturbative, is computed using 
only continuum perturbation theory and has been 
carried out to higher orders, which is well be- 
haved. This method has been shown to be 



quite successful in reproducing results obtained 
by other methods, such as using chiral Ward Iden- 
tities (WI) PJ. The method has also been suc- 
cessfully applied to determine renormalisation co- 
efficients for various operators using the Wilson 
|4l5lo17j and staggered actions [Bj, domain- wall 
fcrmions as well as the quark mass renormal- 
isation constant for overlap fermions ^U]- The 
purpose of the current work is to study the ap- 
plication of this non-perturbative renormalisation 
procedure to the renormalisation of the quark 
field and the flavour non-singlet fermion bilinear 
operators in the case of overlap fermions. 

We shall use Neuberger's overlap fermions [TTj 
which have lattice chiral symmetry at finite cut- 
off. As a result, many chiral-symmetry rela- 

EE 



tions arc still valid [12113) and the quark prop- 
agator preserves the same structure as in the 
continuum. The use of the overlap action entails 
many theoretical advantages [15llrij . such as no 
additive mass renormalisation, no 0(a) error, and 
no mixing among operators of different chirality. 
The latter is very helpful for computing weak ma- 
trix elements. 
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In this section we review the nonperturba- 
tive renormalisation method of Ref. [U, which 
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we use to compute the renormalisation con- 
stants of quark bilinears in this paper. The 
method imposes renormalisation conditions non- 
perturbatively, directly on quark and gluon 
Green's functions, in a fixed gauge, for example 
the Landau gauge. 

Let S (x, 0) denotes the quark propagator on a 
gauge-fixed configuration from a source to all 
space-time points x. The momentum space prop- 
agator is defined as the discrete Fourier transform 
over the sink positions 



S (p, 0) = exp (— ip ■ x) S (x, 0) . 



(1) 



We use periodic boundary conditions in spatial 
directions and an anti-periodic boundary condi- 
tion in time direction. The dimensionless lattice 
momenta are 

2tt N s 2vr 1 N t 

for an N% x Nt lattice, where rij (nt) may in prin- 
ciple lie in the range — > N s (N t ). In practice, 
however, only a subset of this range is used. In 
this paper the momentum range was restricted 
to those momenta for which n, = 0,1,2 and 
nt = 1,2,3,4. 

We also define the square of the absolute mo- 
mentum, 

3 

(pa) 2 = a 2 ^p^p^ . (3) 

2.1. Three Point Function 

Consider the two-fermion operators 

O r (x) = ^(x)T^ , (4) 
where F is the Dirac gamma matrix 



r e {1,7ai>75,7m75,0>} , 



(5) 



and the corresponding notation will be {S, V, P, 
A, T} respectively. The three point function with 
the operator insertion at position and the prop- 
agator from y to and then from to £ is given 
by 

G o (x,0,y) 



= (V(z)O r (0)^(y)) 
= (S(x,0)TS(0,y)) , 



where S(0,y) is the quark propagator from y to 
0. It is the inverse of the Dirac operator. 

The Fourier transform of the three-point func- 
tion is given by 

G ( P a,qa) = J d^xd^ye^'-^Goix^y) 

= <5(p,0)r( 75 5t( (? ,0)75)). (7) 

From this, one can define the amputated three- 
point function 

A (pa, qa) = S(pa)~ 1 G (pa 1 qa)S(qa)^ 1 , (8) 
where 

S(pa) = (S(p,0)) , (9) 

which is translational invariant and a 12 x 12 ma- 
trix. 

Finally, we define a projected vertex function 



T (pa) 



-L— Tr (Ao(pa,pa)P c 



(10) 



Tr(P£, 

where Po = T is the corresponding projection 
operator. 

2.2. RI/MOM Renormalisation Condition 

The renormalised operator 0(fj,) at scale y, is 
related to the bare operator 



0(n) = Z Q (fJ-a, g(a))Q(a) . 



(11) 



and the renormalisation condition is imposed on 
the vertex function To(pa) at a scale p 2 = /i 2 , 



ro,rcn(pa)|p2 =A 



Z ( 



T (pa)\ p 2^2 = 1 , (12) 



to make it agree with the tree-level value of unity. 
Here is the field or wavefunction renormalisa- 
tion 

V'rcn = Z]p 2 lj) . (13) 

In order to obtain the renormalisation constant 
Zo for the operator O, one needs to know Z^ 
first. It has been suggested that Z^ be obtained 
from the renormalisation of vector or axial- vector 
currents. For example, if one uses the conserved 
vector current, then one expects Z v c = 1. There- 
fore, from the renormalisation condition Eq. I|12fl . 
one obtains, 



(0) 



1 



(14) 
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However, in this work, we will obtain di- 
rectly from the quark propagator. It can be de- 
fined from the Ward Identity (WI) as 



-ii-Tr 
12 



dS(pa)- 



(15) 



To avoid derivatives with respect to a discrete 
variable, we have used 



1 T ^J2^i,4l^(P^ a ) s (P a ) 1 



12 



(16) 



which, in the Landau Gauge, differs from Z^, by 
a finite term of order a 2 |17| . and the difference 
is less than 1% at the typical scale p ~ 2 GeV. 
It has been pointed out in Ref. that the sys- 
tematic error due to the definition of lattice mo- 
mentum p is much larger than 1%. We will not 
distinguish Z'^ and Z^ later on. 

3. NUMERICAL RESULTS 

We work on a 16 3 x 32 lattice with lattice 
spacing, a=0.093 fm. The gauge configurations 
are created using a tadpole improved plaque- 
tte plus rectangle (Liischer-Weisz ^B]) gauge ac- 
tion through the pseudo-heat-bath algorithm. A 
total of 50 configurations are considered. The 
lattice parameters are summarised in Table ^ 
The lattice spacing is determined from the static 
quark potential with a string tension ^fa = 
440 MeV 02]. 

Table 1 

Lattice parameters. 



Action 


Volume 


Samp 


(3 a (fm) 


Improved 


16 a x 32 


500 


4.80 0.093 



The gauge field configurations are gauge fixed 
to the Landau gauge using a Conjugate Gradi- 
ent Fourier Acceleration [201 algorithm with an 
accuracy of 9 = J2 \dp,A^(x)\ < 10~ 12 . We use 
an improved gauge-fixing scheme |21| to minimise 
gauge-fixing discretisation errors. 

The massive overlap Dirac operator Jl] is de- 
fined so that at tree-level there is no mass or wave- 
function renormalisation |22j . 

D(mo) =P+^Y + {P- ^)75e(ff) . (17) 



Here e(H) = ff/V H 2 is the matrix sign func- 
tion and H is taken to be the hermitian Wilson- 
Dirac operator, i.e., H = j^Dy?. Here Dw is the 
usual Wilson fermion operator, but with a neg- 
ative mass parameter — p = 1/2/t — 4 in which 
k c < k < 0.25. 

Our numerical calculation begins with an eval- 
uation of the inverse of D(mo) for each gauge con- 
figuration in the ensemble. We use a 14th order 
Zolotarev approximation |23| to the sign matrix 
function e(H w ), and in the selected window of 
x G [0.031, 2.5] of e(x), the approximation is bet- 
ter than 3.3 x 10~ 10 24 . In the calculations, 
k = 0.19163 is used, which gives pa = 1.391. 
We calculate 15 quark masses by using a shifted 
version of a Conjugate Gradient solver. The 
bare quark masses mod are chosen to be 0.0250, 
0.0278, 0.0334, 0.0389, 0.0445, 0.0501, 0.0584, 
0.0668, 0.0834, 0.1001, 0.1252, 0.1669, 0.2086, 
0.2503, 0.2921. The mo values in physical units 
are 53, 59, 71, 83, 94, 106, 124, 142, 177, 212, 
266, 354, 442, 531, and 620 MeV respectively. 

3.1. Axial and vector currents 

Let us consider first the vector and axial vector 
currents. Since each obeys a chiral WI [21] their 
renormalisation constants are finite. In Fig. Q 
we show r,4 and Ty, calculated in the RI scheme 
for different bare quark masses as a function of 
the lattice momentum {pa) 2 . We find they are 
weakly dependent on the mass, and almost scale 
independent for (pa) 2 > 0.5. 

For the overlap fermion, which satisfies the 
Ginsparg- Wilson |26j relation and preserves chi- 
ral symmetry on finite lattice. In Ref. [2U it has 
been proved that the renormalisation constants 
for naive vector and axial vector currents satisfy 



Za = Z\ 



(18) 



In Fig. [21 we show the quantities Y a — and 
and !(r^-|-ry), after extrapolating to the chiral 
limit. We can see from upper part of the figure 
that there is no signal of effects from chiral sym- 
metry breaking, since Y a — Yy is tending to zero 
at moderate and high momenta. The effects of 
spontaneous chiral symmetry breaking are visible 
at low momenta and they are damped at higher 
momenta. At the momenta of interest, there also 



4 



1.3 

1.2 

1.1 

: 1.0 

0.9 

0.8 t- 

0.7 
0.0 

1.4 
1.3 I 
1.3 
■ 1.1 
1.0 
0.9 



||M!;!M!!MI! M 



m a = 0.0334 
m a = 0.0501 

rn a = 0.0B34 
m 'a = 0.1669 
m a = 0.2503 



0.5 



1 .0 



(pa) 2 



1.5 



2.0 



2.5 



0.8 



!' 



m a = 0.0334 
m q a = 0.0501 

m q a = 0.0834 
m q a = 0.1669 
m q a = 0.2503 



0.0 



0.5 



1 .0 



(pa) 2 



1.5 



2.0 



2.5 



Figure 1. The projected vertex function F defined 
in Eq. (1101 for vector and axial vector currents with 
different bare quark masses. By the renormalisation 
condition (Eq.(^J), they are the ratio of renormali- 
sation constant Z-^/Za and Z^jZy respectively. 

seems to be no significant splitting due to non- 
perturbative effects with the difference between 
and IV being less than 1% at (ap) 2 = 0.7 and 
smaller for momenta above this. In the lower part 
of Fig. 13 we plot \(Ta + Ty) against (ap) 2 , and 
it can be used for the extraction of both Za/Z^ 
and Zy 1 1 Z^p to increase the statistical accuracy. 

3.2. Pseudoscalar and scalar densities 

The pseudoscalar and scalar densities differ 
from the axial and vector currents in that their 
renormalisation is not scale independent. For 
overlap fermions with chiral symmetry the RI 
scheme preserves the well known MS relations 



and Z n 



1 

Zs 



(19) 



and the quantities Z$/Zp, ZgZ rn and Z m Zp are 
expected to be scale independent. 
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Figure 2. Checking Za = Zv, the upper plot is 
(Fa — IV) versus (pa) 2 , where we linear extrapolate 
to the chiral limit (mo = 0). The lower plot is i 
(Fa + IV) versus (pa) 2 . The relation Za = Zv is 
valid at moderate to large (pa) 2 . 

Fig. Eland Fig. ^show Ts and Tp versus (pa) 2 
with different masses. Unlike in the case of the 
vector and axial vector currents, they are strongly 
dependent on the quark mass tuq for light quarks. 
When we want to extract the chiral limit value, 
we cannot use a simple linear or quadratic fit. In 
Fig.|3|and Fig. 01 we can clearly see the pole effect 
in the small (pa) 2 region. 

As discussed in detail in Ref. 0, due to the 
contribution from the zero-modes, we may fit the 
Ts to the form 



r, = 



Cl,S 

(am )' 



- C2,s + C3,s(am ) 5 



(20) 



at a fixed momentum with Z^jZg being given by 
C2,s- The fitted value of C2,s is plotted in Fig.Elas 
Z^jZs- We find that the coefficient of the pole 
term, ci } s is very small, between 10~ 3 to 10~ 4 , 
this means that the pole term only provides a 
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Figure 3. Fs versus (pa) 2 with different masses. Here 
one can see the strong mass dependence and non- 
monotonic behaviour. 
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Figure 4. Fp versus (pa) 2 with different masses. 
Here also we have strong mass dependence and pole 
behaviour. 

large contribution at very small mo, i.e., near the 
chiral limit. 

For Tp, we fit to the form 



r 



Cl,P 



P.latt 



C2,P 



(am ) 2 (am ) 



+c 3i p+c^ P (am ) 2 , (21) 



with C3,p being equal to Z^fZp. The quadratic 
mass pole is due to zero-mode effects in (qq) . The 
fitted value of C3,p is plotted in Fig.0 clS Zip /Z P . 

As the evidence for the above fitting form, 
we consider the resulting values for Z^/Z$ and 
Z^/Zp. Fig. [3] shows a comparison between the 
extracted values of these two quantities. As chiral 
symmetry would predict for Z^/Zg and Z^jZp, 
the two quantities coincide at moderate and large 
momenta. This provides an excellent test of both 
the fitting method to extract the poles and the 
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Figure 5. Checking of Zp = Zs after mass pole 
subtraction. 
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Figure 6. FV versus (pa) 2 with different masses. It 
shows no significant mass dependence for moderate 
to large (pa) 2 . 

chiral properties of overlap fermions. 

3.3. The Tensor current 

In Fig.|S|we show Tj- versus (pa) 2 with different 
masses. We can see that at moderate and large 
(pa) 2 , Tt is not sensitive to the quark masses. 
The chiral limit value is obtained by a linear fit 
in quark mass as for the cases of vector and axial 
vector currents, and the results will be presented 
in the next subsection. 

3.4. Running of the renormalisation con- 
stants 

The renormalised operators are defined as 



ZoOba 



Or. 



(22) 



Requiring the bare operator to be independent of 
the renormalisation scale /j, 2 gives the RG equa- 
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Figure 7. versus (pa) 2 in the chiral limit before 
(labelled as "bare") and after (labelled as "SI") three 
loop perturbative running, such that they coincide at 
(pa) 2 = 1. The latter is almost (pa) 2 independent 
after (pa) 2 > 0.6. The slope versus (pa) 2 is about 
0.01 and is interpreted as an 0(a 2 ) effect. 
1.1 



1.6 



1.2 



1.0 



0.E 





i 


1 




D bare 






* SI 




a D a □ s □ 

.:hi 


x s « a x * * 






Ho 


H □ Q 








E ^_ 







0.0 



0.5 



1.0 1.5 
(pa) 2 



2.0 



2.5 



Figure 8. The same as Fig. ||2J. Here is Fs versus 
(pa) 2 . The slope of SI versus (pa) 2 is about 0.02. 

tion, 



1 2 <^0 



Ore 



^Orcn. (23) 



d M 2 Z(j 2 

The solution of Zo(n 2 ) can be written in the fol- 
lowing form 



Zo (M 2 ) 



Co (n 2 ) 
Co (m' 2 ) 



(24) 



The coefficient Co (a* 2 ) can be found in Appendix 
B of Ref. |Hj and references within. 

In this work, the value of a s was calculated at 
three loops using a lattice value of Aqcd taken 
from Ref. 28, as A QC d = 238 ± 19MeV . 



Figure 9. Ft versus (pa) 2 in the chiral limit before 
and after two loop perturbative running. The latter is 
almost (pa) 2 independent after (pa) 2 > 0.6, the slope 
versus (pa) 2 is about 0.01. 

Both Za and Zv should be scale independent, 
but this is not the case for Z$. Fig. shows 
both r^4 and the scale invariant (SI) quantity (the 
data after removing the renormalisation group 
running) calculated as described above: 



r s / ((apf)=r A ((ap) 2 )/a k ((apf 



(25) 



The quantity C$ is normalised so that C^(l) = 1. 
As can be seen, in this case the renormalisation 
group running is coming from Z^ alone. It is 
small, but it actually improves the scale indepen- 
dence of the data. The remain scale dependence 
of this data is very small and a plausible explana- 
tion for this is a small (ap) 2 error. Indeed, when a 
linear fit of the SI data versus (ap) 2 is performed, 
for 0.8 < (ap) 2 < 2.0, the gradient is w 0.01. 

In the case of Z^/Zs, the two renormalisation 
constants are both running with (ap) 2 . Fig. [S] 
shows both Ts and the scale invariant (SI) quan- 
tity after three loop running. The modification to 
the raw data, which is taken after the mass-pole 
has been subtracted gives a satisfactory result. 
The linear fit of the SI data versus (ap) 2 in the 
range of 0.8 < (ap) 2 < 2.3, gives a gradient of 
about 0.02. 

The SI result of Tt which comes from the two 
loop running of Zt |29| and three loop running 
of Z^i, is plotted in Fig. [5J The "bare" data is 
taken by a simple linear extrapolation to the chi- 
ral limit. The linear fit to the SI data gives gra- 
dient about 0.01. 
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Figure 10. A(p) in the chiral limit and its scale 
invariant counterpart versus (pa) 2 as a + b (pa) 2 . 
The fitting range is [0.70, 2.32], which gives Z$ about 
1.133 ±0.006. 

With the interpretation that the remaining 
scale dependence is due to 0(ap) 2 effects, the 
correct way to extract the renormalisation coeffi- 
cients is to first construct the SI quantity as de- 
scribed above, and then fit any remaining scale 
dependence 0] to the form 



V = ci + c 2 (ap) 2 , 



(26) 



for a range of momenta that is chosen to be 
"above" the region for which condensate effects 
are important. We choose the range 0.8 < 
(ap) 2 < 2.3. We then apply the renormalisation 
group running formula (an inverse operation of 
Eq. ) to ci in Eq. J2BJ) to get the ratio of 
the renormalisation constant which is now free of 
(ap) 2 error in RI scheme. 

3.5. Extracting from the propagator 

Here we will use Eq. (|16|) to extract the wave 
function renormalisation constant Z^. We will 
use two different definitions of lattice momenta, 
i.e. , the discrete lattice momentum defined in 
Eq. (J2J and the kinematic lattice momentum 
which was introduced in Eq. (A9) of Ref. |3()| . 

3.5.1. Discrete lattice momentum p^ 

We here use the discrete lattice momentum p^ 
in Eq. ifTBT) to calculate Z^. By using Eq. ijToj) . we 
can get the corresponding quantity of ZL (here we 
use the symbol A(p) for ZL). The resulting A(p) 
is strongly dependent on the scale. We use three 
loop pcrturbative running in RI scheme and fi- 
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Figure 11. Z^ calculated from Eq. 1151 by using 
the kinematic lattice momentum g M instead of the 
discrete lattice momentum p^. 

nally extract the value of Z^, from A(p) by fitting 
A(p) to Eq. (12(it for a range of momenta for which 
the condensate effects are unimportant. Here we 
chose the range (ap) 2 C (0.70,2.32), which cor- 
responds to p 6 (1.78, 3.23) GeV. ci can be iden- 
tified as the SI value of Z^. In Fig. ^|we show a 
linear fit for A(p) versus (pa) 2 as in Eq. H2ri|) . 

3.5.2. Kinematic lattice momentum 

Another way to calculate Z^ is to use the kine- 
matic lattice momentum q^ defined in Eq. (A9) 
of Ref. [3U]. The calculations are similar to those 
using p tl . The result is plotted in Fig. ^2 The 
fitted value of Z^ is 1.179(5), with the difference 
to the value calculated by p M being less than 4%. 
In both calculations, we attribute the remaining 
(pa) 2 dependence to the possible 0((pa) 2 ) error 
and use the linear fit Ea.f l26|l to remove it. When 
using the kinematic lattice momentum g M to cal- 
culate Z^, the 0((pa) 2 ) errors are much smaller 
than in the case of using the discrete lattice mo- 
mentum p^. 

3.6. Matching to MS scheme 

In order to confront experiment, one frequently 
prefers to quote the final results in the MS 
scheme at a certain scale. For light hadrons, the 
popular scale is 2 GeV. 

The perturbative expansion of the ratio 
Z MS /Z RI to two-loop order is given by 0] 



7MS 



R 



§r = 1+ g (z » ) (.) +( g )2(z 







.(27) 
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The numerical values of the matching coefficients, 
Z ( Q 1] and Z ( 2) in Eq. |23) for Z$ and Z s can be 
found in Appendix C of Ref. we have not 
found the numerical value for Zt in the literature. 
The final results for the renormalisation constants 
are listed in Tabled 

Table 2 



Final Z-factor results 



Z — facte 


ir 


RI (2 GeV) 


MS at 2 GeV 


Za/Z^j 
Zs/Z^j) 
Zt/Z^ 
Z^(from 
Zip(from 


P) 

?) 


0.924±0.004 
0.739±0.003 
1.009±0.002 
1.134±0.006 
1.180±0.005 


0.928±0.004 
0.842±0.004 
1.012±0.002 
1.130±0.006 
1.175±0.005 



4. SUMMARY 

In this work we performed the non-perturbative 
renormalisation of composite operators with over- 
lap fermions. The main results are: in MS 
scheme at \i = 2 GeV, Z^, = 1.153(5)(22), Z s 
= 0.970(4)(20), Z A = 1.070(5)(21), and Z T = 
1.167(3) (23), where the first error is statistical 
and the second error comes from using the two 
different lattice momenta to get Z^. 

Support for this research from the Australian 
Research Council is gratefully acknowledged. 
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